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Abstract. The unique core- mass - radius relation for gi- 
ants with degenerate helium cores enables us to recon- 



Models for the formation of close double white dwarfs 
envision two standard scenarios to produce these systems 



struct the evolution of three observed double helium white ( 


Futukov & Yungelson 1981 




[ben & Tutukov 1984 




Web- 


dwarfs 


with known masses of both components. bink 1984 




Futukov & Yungelson 1988 




Han et al. 1995; 


The last mass transfer phase in their evolution must [ben ct al. 1997 




Ban 1998). In the first 


scenario, two low- 



have been a spiral-in. In the formalism proposed by Web 



bink (1984), we can constrain the efficiency of the deposi- 
tion of orbital energy into the envelope to be 1 <, a 55 6, 
for an envelope structure parameter A = 0.5. We find that 
the two standard mass transfer types (stable mass transfer 
and spiral-in) are both unable to explain the first phase 
of mass transfer for these three binaries. 

We use a parametric approach to describe mass trans- 
fer in low-mass binaries, where both stars are of compara- 
ble mass and find that the orbital characteristics of the ob- 
served double helium white dwarfs can be well reproduced 
if the envelope of the primary is lost with ~1.5 times the 
specific angular momentum of the initial binary. In this 
case no substantial spiral-in occurs. 

Key words: stars: white dwarfs - stars: mass loss - bi- 
naries: close - stars: evolution 



1. Introduction 

The long lasting problem that we observe many double 
stars which are expected to form close pairs of white 
dwarfs, but yet that of the observed white dwarfs not one 
seemed to have a close white dwarf companion, was solved 
by the discovery of such pairs, starting wit h L870-2 (= 
WD 0135+052) in 1988 QSaffer et al. 1988| ). In total 14 
close detached binary white dwarfs are known at present, 
see Table |l[ The fact that six of these systems have their 
orbital period and the masses of both components deter- 
mined provides an opportunity to test binary evolution 
theory in detail. 



mass (M SS 2.3 Mq) stars evolve through two stages of 
spiral-in. The first spiral-in will shrink the orbit, so the sec- 
ond spiral-in happens in a binary with a much smaller or- 
bital separation. The Roche lobe filling giant (secondary) 
now has a small radius and therefore a small core. The 
white dwarf that is formed last is thus less massive than 
its companion, with a mass ratio mbright/^dim Si 0.55 (see 
e.g. [Sandquist ct al. 200C ). 

In the second scenario, the first phase of mass trans- 
fer is stable; the second phase of mass transfer is again 
a spiral-in. The evolution of the orbit and the growth 
of the core during the first, slow phase of mass trans- 
fer depend on the amount of mass and angular momen- 
tum that is lost from the system. If the evolution in 
this phase is conservative, t he expected final mass ratio 
Wbright [rndim ~ 1-14 — 1.18 ( Tutukov fc Yungelson 1988 ; 
fcandquist et al. 2000|) . 



All white dwarfs in close pairs known today have low 
masses (M <, 0.5Mq). Note, however, that the inaccuracy 
of the mass determinations is as large as ~ 0.05Mq due 
to uncertainties in model atmospheres and cooling curves 
for white dwarfs (see e.g. Napiwotzki et al. 1999 ). These 
low masses suggest they are helium white dwarfs, but it 
cannot be excluded a priori that white dwarfs with masses 
<; 0.35 Mq are so called hybrid white dwarfs, i.e. having 
small CO cores and relatively th ick (~ 0.1 Mq) helium 
envelopes ( [ben fc Tutukov 1985| ). For the most massive 
ones (M ^ 0.45 Mq), this is even the only option, since 
helium white dwarfs must have a mass below 0.46 Mq 
(Sweigart et al. 199C). For the less massive ones the prob- 
ability to form hybrid white dwarfs is 4 - 5 times lower 



than to form helium white dwarfs (Nelemans et al. 2000) 
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Table 1. Parameters of known close double white dwarfs 
with towd denoting the mass of the brightest white dwarf 
and A/\vd denoting its companion. For references see 
Maxted fc Marsh (1999I ) 1 and |Maxted et al. (200C| ). 



WD 


P(d) 


777 wn /iVT^-\ 
' ' " W i-J / "CO 


Mwn /Mrs 




1 ^fi 






0136+768 


1.407 


0.44 


0.34 


0957-666 


0.061 


0.37 


0.32 


1101+364 


0.145 


0.31 


0.36 


1204+450 


1.603 


0.51 


0.51 


1704+481 


0.145 


0.39 


0.56 


1022+050 


1.157 


0.35 




1202+608 


1.493 


0.40 




1241-010 


3.347 


0.31 




1317+453 


4.872 


0.33 




1713+332 


1.123 


0.38 




1824+040 


6.266 


0.39 




2032+188 


5.084 


0.36 




2331+290 


0.167 


0.39 





1 F or WD 0136+768 we give the masses of com ponents af- 
ter Bergeron et al. (1992 ), correcting a misprint in 



Marsh 



199£ ) 



Maxted 



2. Reconstructing the binary evolution 

Because of the unique core-mass - radius relation for gi- 
ants with degenerate helium cores (Refsdal & Weigert 
197Cl|)7|we can reconstruct the mass transfer phases in 
which helium white dwarfs are formed. The mass of the 
brightest star in WD 1101+364 (0.31 M©) indicates that 
it is a helium white dwarf. In WD 136+768 and WD 0957- 
666 it cannot be excluded that the brightest stars are hy- 
brid white dwarfs. The low mass of the dimmer compan- 
ions in these three systems indicates that those are all 
helium white dwarfs. The white dwarfs in WD 0135-052 
and WD 1204+450 are formally inconsistent with being 
helium white dwarfs, so we will not include them in the 
discussion anymore. WD 1704+481 probably consists of a 
helium white dwarf and a dimmer CO white dwarf, but 
because CO white dwarfs cool faster than helium white 



dwarfs ( Driebe et al. 1998 ), it is not clear which of the 
white dwarfs was formed most recent so we cannot use 
this system in our present study. 

Assuming that the mass of the white dwarf is equal 
to the mass of the core of the giant at the onset of the 
mass transfer, the radius of the progenitor of a helium 
white dwarf can be calculated from the core- mass - radius 



relation given by [ben & Tutukov (1985) 
R w 10 3 - 5 M 4 



(1) 



(R and M c in solar units). This equation is in good agree- 
ment with other equations describing this dependence for 
giants (e.g. Webbink et al. 1983| ). The mass of the white 
dwarf progenitor is however not known, since the above 
relation is independent of the total mass of the star. How- 
ever the mass of the giant must be in the range 0.8 - 2.3 
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Fig. 1. Periods before the spiral- in phase in which the 
younger white dwarf was formed as function of the mass 
of its giant progenitor. Lines from bottom to top: WD 
1101+346, 0957-666 and 0136+768 



Mq. For less massive stars the main-sequence life time is 
larger than the age of the Galaxy. More massive stars do 
not form degenerate helium cores. 

For the remainder of this article we use the following 
notation: Mi and AfyyD indicate the initial mass of the 
original primary and the mass of the white dwarf that it 
forms, m\ and kiwd represent the same for the original 
secondary. If the secondary accretes mass during the first 
phase of mass transfer, we represent its new mass with 
m(. For the radii of the stars when they become giants 
we use R g and r g for the original primary and secondary 
respectively. With P U P 1U and P we indicate the initial 
period, the period after the first phase of mass transfer, 
and the current period of the binary. 

3. Last mass transfer: spiral-in 

Using Eq. ([!]) we calculate the radii of the progenitors 
of the brightest white dwarfs for the three double helium 
white dwarfs. Since we know the mass of the white dwarf 
that orbited this giant and may reasonably assume that 
it did not accrete anything during the spiral-in phase, we 
can calculate the orbital separation at the onset of the 
spiral-in, as function of the mass m[ of the giant, 



>K) = 



r g (m WD ) 
r L (M WD /mO' 



(2) 



where r^ = R^/a is the dimcnsionlcss Roche lobe radius, 
given e.g. by Eggleton (1983) and we assume r g = i?L- 
This is shown in Fig. |1|, where we use Keplers 3 rd law 
to compute the period from the orbital separation. Com- 
paring the periods in Fig. |l| with the observed periods in 
Table [j], we see that in the last mass transfer phase the 
orbital separation must have reduced dramatically. This 
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Fig. 2. The parameter aX for WD 0957-666, 1101+364 
and 0136+768 assuming that the brightest component is 
a helium white dwarf and their orbital periods did not 
change since the end of the spiral-in stage. Lines from bot- 
tom to top are for WD 0957-666, 0136+768 and 1101+364. 

can only be accomplished if the last mass transfer was a 
spiral-in. 

In a spiral-in, the envelope of the giant is expelled at 
the expense of the orbital energy of the binary. Balancing 
the binding energy of the envelope of the giant with the 
difference in orbital energy (Webbink 1984) one finds 



Mwd (m| - rawp) 
A r„ 



M W D mwD Mwd m[ 



2 of 



2 a„ 



(3) 



The parameter A depends on the structure of red giant 
envelo pe. T he usual assumption is that A = 0.5 (Dc Kool 



et al. 1 987 ). The parameter a represents the efficiency of 
the deposition of orbital energy into the common envelope. 
To reduce the number of parameters, the product a A is 
treated as a single parameter in the remainder of this ar- 
ticle, but it should be noted that both A and a will in 



reality Ibe functions ot the evolutionary stage of the stars. 

Applying Eq! (||) and Fq. (||), we find aA as a lunc- 
tion of m[ for the three systems considered. We plot this in 
Fig. ||, where we assume that the current periods are equal 
to the post spiral-in periods. This may not be the case in 
general since close orbits like the those of WD 0957-666 
and WD 1101+364 will decay due to the loss of angular 
momentum by gravitational radiation. However, the esti- 
mated ages for these white dwarfs (~ 10 7 yr (Moran ct al. 
199'iirknd ~ 10 9 yr (using the cooling curves of pricboj 



et al. (1998)) respectively) are short compared to the or- 
bital decay time scale. 

For the remaining white dwarf pairs listed in Tabic [l] 
the mass of only one component is known. We assume here 
that it is the last formed component we observe. Low-mass 
white dwarfs may have thick hydrogen envelop es which 
make them cool very slowly ( Driebe et al. 1998 ) and the 
situation in which the older white dwarf is really observed 
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Fig. 3. Periods before the spiral-in for all double white 
dwarfs in which the last formed is a helium white dwarf. 
Limits allow for the (unknown) mass of the progenitor of 
this white dwarf and the mass of the unseen companion. 
Solid lines are for a companion mass between 0.2 and 0.65 
Mq. Dotted line gives the limit for a companion mass of 
0.65 - 1.4 M 



a priori. 



However, as we show in an 



cannot be excluded 
forthcoming paper (Nelcmans et al. 2000), in the majority 



of binary white dwarfs we indeed observe the youngest of 
the two dwarfs. 

From Figs, [l] and || we see that we find a range of 
P m 's and aA's for WD 0957-666, 1101+364 and 0136+768 
where the mass of the second white dwarf is known. For 
the remaining systems we can also compute a range of 
P m 's and aA's by determining the ranges for all possible 
masses of the unseen companion. We do know that the 
companion almost certainly is another white dwarf, so the 
mass of this object must be betwee n 0.2 and 1.4 Mq and 
most probably even below 0.65 Mq (Iben et al. 1997; Han 



1998; Nelemans et al. 2000). Having also in mind that the 
mass of white dwarf progenitor is in the 0.8 - 2.3 Mq 
range we can derive the possible ranges of intermediate 
periods and aA's for all double white dwarfs, as is shown 
in Figs. | and g 

Two conclusions can be drawn from Figs, [l] to ||: 

1. The efficiency of the energy deposition into the com- 
mon envelopes a must be high. From model calcula- 
tions of stellar structure we know that A ~ 0.5 — 1.0. If 
we assume that a does not depend on the evolutionary 
state of the giant or the combination of the masses of 
the giant and the white dwarf, the parameters of all 
double white dwarfs with two observed helium com- 
ponents can be reproduced with the same a ~ 4. The 
only exception is WD 0957-666 for w hich the efficiency 
appears to be much lower (see Sect. |6.2| for a different 
solution) . Note that an error in the masses of the white 
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Periods after the first mass transfer 
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Fig. 4. The parameter aA for the same cases as in Fig. || 

dwarfs of 0.05 Mq translates to an error in the value of 
a A of a factor ~1.5. However even this does not bring 
the value of a to 4 for WD 0957-666. 
Since Eq. (||) only considers a rough energy budget, the 
conclusion a > 1 could simply mean that we do not 
calculate the energy accurately. It does mean that the 
orbital energy deposition into common envelope has 
to be highly efficient. It could also mean that sources 
other than the orbital energy contribute to the process 
of common envelope expulsion ( Iben fe Livio 1993| ). 
E.g. it is possible that the envelope is partially removed 
before Roche lobe contact by an enhanced stellar wind 
due to tidal interaction between the giant and the com- 
panion ( |Tout fc Egglcton I988| ) , yielding a lower value 
of aX. 

2. The immediate progenitors of the known close double 
white dwarfs (i.e. the white dwarf + giant binaries) 
all had rather wide orbits (25 - >500 days). This has 
important consequences for the understanding of the 
first phase of mass transfer. 



4. The first mass transfer 

We compute the evolution of the binary parameters in 
the first mass transfer, where we start from the initial 
binary and evolve it forward according to the two stan- 
dard scenario's. The resulting periods should be equal to 
the intermediate periods we reconstructed in the previous 
section. 



4.1. Spiral-in 

In the case when the first mass transfer was also a spiral-in, 
we can compute the period after the spiral-in by making 
a (standard) assumption that aA is constant (i.e. here, as 
derived above aA = 2). The initial separation at which the 
primary fills its Roche lobe is again determined by Eq. ([]]) 
and a value for Mj between 0.8 and 2.3 Mq. Applying 
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Fig. 5. Periods after the first phase of mass transfer (P m ) 
as function of the mass of secondary at this time (m(). 
Top three lines are periods needed to explain the mass of 
the second formed white dwarf (see Fig. Q). Middle three 
lines give the maximum period would the first phase of 
mass transfer be an Algol phase, lower three ones are for 
the case the first mass transfer phase was a spiral-in. Solid 
lines for WD 0957-666, dashed lines for WD 1101+346 and 
dotted lines for WD 0136+768. 



Eq. (|3|), we find a m (and thus P m ) as a function of Mj 
and to;. The maximum period after the spiral- in is given 
for the case Mj = mi, which we plot in Fig. || (bottom 
three lines). These periods are clearly much smaller than 
the periods derived in the previous section (see Fig. [l]) . We 
reproduce the latter on a logarith mic scale as the top three 
lines in Fig. |[ In a recent article Sandquist et al. (2000 ) 
came to the same conclusion, based on the observed mass 
ratios and also concluded that it is very hard for bina- 
ries containing two low-mass stars to survive two spiral-in 
phases. 

It may be argued that aA is different in the first mass 
transfer phases, because the companion now is a main- 
sequence star instead of a white dwarf. We can, just as in 
the case of the second mass transfer determine the value of 
aA that is required to get the right period after the mass 
transfer. If we do this, however, we find — 15 ^ aA < —5 
clearly out of the allowed range. This means that a spiral- 
in phase in the first mass transfer phase for these systems 
is ruled out. 



4-2. Stable mass transfer 

For the alternative scenario, to start with a short period 
zero-age binary and evolve through a phase of stable mass 
transfer, the zero-age orbital period should be in a narrow 
interval, such that the primary fills its Roche lobe at the 
moment it still has a radiative or at least shallow convec- 
tive envelope. For these stars mass transfer in which the 
donor stays in thermal equilibrium is possible if the ini- 



5 



tial mass ratio q\ = m ; /Mi > 0.8 ( [Tutukov et al. 1982| ). 
For more extreme mass ratio's; q\ ^ 0.3, mass transfer 
can proceed on a timescale in between the thermal and 
the dynamical timescale of the donor (e.g. trial computa- 
tions for a system with Mi = 2.3MQ,mi = 0.8Mq and 
di = 12.8i? show that log Mi <> -5.5M yr" 1 ; A. Fe- 
dorova, private communication). However for mass ratio's 
q\ 0.5, the thermal timescale of the accretor is much 
longer than the mass transfer timescale and the accretor 
is expected to expand, leading to a contact configuration 



Table 2. Parameters of wide binaries with giant and white 
dwarf components. Masses in solar units. 



Name 



P (d) M m W D 



ref 



S1040 42.; 
AY Cet 56.; 



L5 022 [Landsman et al. ( 19971) 



2.2 



0.25 



Simon et al. (1985) 



which s unstable (e.g. Kippcnhahn & Meycr-Hofmeister 



1977|)7rnlcss the secondary is a convective low-mass main- 
sequence star (Webbink 1977). 

Since we want to explore the limits of the Algol sce- 
nario, we allow all values of qi > 0.3 and calculate the 
possible periods after the Algol phase. 

We assume the mass transfer is conservative and for 
a given total mass of the system M to t = M w< \ + m[, the 
initial masses are M; = M t ot/(l + qC) and to; = ftMi. 
We consider systems in which the primary fills its Roche 
lobe just before the star develops a relatively deep convec- 
tive envelope (the outer 50% of the mass is convective). 
The radii of such stars are ob tained by a fit to stellar 
models by pVlengel et al. (1979| ): i? max w 2.4 M; 156 (solar 
uni ts). Using the ap proximation to the Roche lobe given 
by Paczynski (1967| ) and the equation for the change in 
period for conservative mass transfer we find the period 
after the Algol phase is 



1.38 



M, 



7.81 



Qi 



(M wd m[) 3 (l + qi ) 



7.81 



days. 



(4) 



This equation has a maximum for qi = 0.62. We com- 
pute these maximal periods as function of m[, where 
of course Mi and m; are limited by Mi > 0.8Mq, be- 
cause the primary has to evolve off the main-sequence and 
m[ < 2.3Mq, because the secondary must still develop a 
degenerate helium core. 

The resulting periods for q = 0.62 are shown as the 
middle set of lines in Fig. 0, where we actually used the 
equation for the Roche lobe given by [Eggleton (1983] ), 
which is better in this mass ratio regime. The obtained 
periods are clearly not long enough to explain of the origin 
of WD 0957-666, 1101+364 and 0136+768. 

One could argue that the assumption of conservative 



evolution is not correct. Han (1998), for example, in his 
'best' model assumes 50% of the mass is lost with the 
specific angular momentum of the donor. This could yield 
wider orbits than obtained above. Looking at his Fig. 7 
one sees that he can indeed explain WD 1101+364, but 
not WD 0957-666 and especially not WD 0136+768, be- 
cause for these systems the masses of the last formed white 
dwarfs (i.e. the periods after the Algol phase) are too 
large. To get higher masses one needs to lose less angular 
momentum with the mass, leading to even wider orbits. 
However, analysis of observed low-mass Algols, which are 
binaries currently in this stage of stable mass exchange 



e.g. Refsdal et al. 1974 ; Massevich fc Yungelson 1975 



|Giannuzzi 1981 ; Iben fe Tutukov 1984 ; Kraicheva et al 



1986| ; [Maxted fc Hilditch 1996| ) shows that their periods 
are smaller than would be expected with Eq. (|J) . This sug- 
gests that descendants of Algol-type systems have orbital 
separations even smaller than in the case of conservative 
evolution, which we assumed for Fig. ^[ 

We conclude that with the above assumptions Al- 
gol evolution in the first mass transfer stage for WD 
0136+ 768, 0957-66 and 1101+364 is ruled out also. 

There are some other observed systems with white 
dwarfs which probably could not be formed through Algol- 
type evolution (Table ||), since white dwarf + giant bina- 
ries with Mwd 0.2 5M(i) end their Algol phase with 
periods below ^30 d. ( Kraicheva et al. 1986| ). In another 
system, HD 185510, with P or b=20.7 day, a giant has a hot 
companion which is classified as a 0.3 Mq sdB st ar from 
its temperature and gravity ( Frasca et al. 19*9^ ). How- 



ever, it also fits the range of temperatures and gravities 
for - 0.25Mq helium white dwarfs ( priebe et al. 1998[ ) 
and the system may be a viable Algol descendant. If it 
is an sdB star, which a priori is less likely, it actually 
matches a scenario similar to one for WD 0957-666 rather 



well (see Sect. QJ2 



5. Unstable mass transfer revised 

Since both standard scenarios for the first phase of mass 
transfer appear to be ruled out, the situation apparently 



is more complex. Tout & Eggleton (1988) (see also Han 



1998) assume that due to tidal effects of the companion, 



the star can lose up to 150 times more mass in the stellar 
wind than without companion. This has of course the de- 
sired effect that the orbit widens before the mass transfer 
starts, and that there is less envelope mass left to be ex- 
pelled, leading to a less dramatic spiral-in. However, if we 
recompute the lines for a spiral-in in Fig. |5| with reduced 
envelopes such that they overlap with the reconstructed 
lines, we find that for WD 1101+364 we need to reduce 
the envelope by 70% and for the others by even more than 
90%. 

In searching for a different solution we start by noting 
that the original spiral-in picture ( [Paczyhski 1976| ) con- 
siders a companion which orbits inside the envelope of the 
giant. The companion experiences drag forces while mov- 
ing in the envelope and frictional effects brake the com- 
panion. In this process orbital energy is transformed into 
heat and motion of the gas and finally into kinetic energy 
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that causes the envelope of the giant to be expelled. This 
picture is very much based on the situation where there 
is a tidal instability which causes the decay of the orbit of 
the companion in systems with a high mass ratios of the 
components (Darwin 1908; Counsclman 1973] ). In the case 
that the common envelope is caused by a runaway mass 
transfer, the common envelope will not look much like the 
equilibrium envelope of the star and worse, the angular 
momentum of the orbit is so large that the common enve- 
lope is in principle easily brought into co-rotation with the 
orbit. At that moment there are no drag forces anymore. 

Since for stars with deep convective envelopes mass 
loss on a dynamical time scale seems, in the current state 
of the art in stellar evolution modelling, inevitable, we 
have to assume that in the progenitors of the observed 
double white dwarfs some kind of common envelope en- 
gulfing the whole system forms. The parameters of the 
observed close binary white dwarfs suggest that this en- 
velope is subsequently lost without much spiral-in. The 
energy to expel the envelope may be supplied by the lu- 
minosity of the giant or by tidal heating, or a combination 
of both. In absence of a detailed physical description we 
will describe the effects of this mechanism in terms of the 
angular momentum balance. 

We compare the pre- and post-mass transfer binaries, 
under the assumption that the envelope of the giant is lost 
completely from the binary (i.e. m[ — m;), and that this 
mass loss reduces the angular momentum of the system 
in a linear way, as first suggested for the general case of 



non-conservative mass transfer by Paczyhski & Ziolkowski 
(19670 



AM 



Aft, 



(5) 



where J\ is the angular momentum of the pre-mass trans- 
fer binary and M tot is the total mass of the binary. The 
change of the orbital period as function of the initial and 
final masses of the components then becomes 



0136+768 



1 101+364 



0957-666 



Pm = { M WD m[ 
Pi ~ I Mmii 



M 



WD 



Mj + m ; 



1-7 



AM 



.(6) 



We can estimate 7 for the three double helium white 
dwarf systems, just as for the case of aX. It turns out that 
all three systems are consistent with a value of 7 between 
~1.4 and ~1.7 (Fig. §). 

Thus the mass transfer from a giant to a main-sequence 
star may in general cither be stable (in the case where 
the giant still has a radiative, or at least not too deep 
convective envelope), unstable, leading to a spiral- in, or a 
process in which the envelope is lost without much of a 
spiral-in. Which systems do and which do not experience 
a spiral-in is related to the mass ratio of the components. 
As can be seen from Eq. for systems with 7 AM w 
Mi (1 + q) the periods already become very small and the 
effect is essentially the same as in the case of a spiral-in. 



1.5 



2.5 



Fig. 6. Derived range of possible values of 7, for the three 
double helium white dwarfs. 



6. Formation of observed systems 

From Figs. | and | it's clear that WD 0136+768, WD 
0957-666 and WD 1101+364 could be formed with aX 
below 0.8. However from Fig. ^ it looks like the extremely 
short-period system WD 0957-66 falls out of the sample 
in the sense that the other systems, as well as systems 
with unobserved companions, are compatible with a value 
of aX ~ 2. 



6.1. Formation of helium white dwarf pairs 

In Fig. ^we show evolutionary scenarios for WD 0136+768 
(left) and WD 1101+364 (middle) in which they consist of 
two helium white dwarfs. We included the effect of stellar 
winds which was not taken into account in the preceding 
discussion. Therefore we use slightly different values for 
7 (1.75 and 1.85 respectively) than the values derived in 
Fig. ^. The difference in the initial mass ratio of the com- 
ponents results in dramatically different orbital periods 
after the envelope of the primary is lost: almost tripling 
in the first case and decrease by half in the second. As 
a result, after the second mass transfer episode the sec- 
ond white dwarf is the more massive member of the pair 
in WD 0136+768 but the less massive in WD 1101+364. 
The large difference in periods and especially masses after 
the first phase of mass transfer, results in rather different 
final orbital periods. 

For WD 0957-666 a scenario in which both components 
are helium white dwarfs can also be constructed but then 
one has to assume that aX is atypically small. However 
we suggest a different scenario. 

6.2. An alternative scenario for WD 0957-666: 
carbon- oxygen white dwarf with helium companion 

There is another solution which allows us to explain the 
origin of WD 0957-666 using aX ~ 2 and 7 « 1.75, like 
for the other two systems. 
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WD 0136+768 




M 2 (M ) P b (d) 

2.12 64,0 

2.12 68.9 
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M,(M S ) 



WD 1101+364 
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WD 0957-666 




M 2 (M ) P b (d) 
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0.31 0.1 
0.31 0.1 
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Fig. 7. Evolutionary scenarios for the formation of WD 0136+768, WD 1101+364 and WD 0957-666 (left to right). 
In all scenario's the primaries lose their envelope after filling their Roche lobe, causing a change in the orbital period 
described by Eq. (|). Whether the orbit widens (WD 0136+768 and WD 0957-666) or shrinks (WD 1101+364) depends 
on the mass ratio. The second mass transfer always results in a spiral-in. For WD 0957-666 we present a scenario in 
which in the first mass transfer a helium star (sdB star) is formed which becomes a hybrid white dwarf only after the 
companion has become a helium white dwarf in the second phase of mass transfer. 



The mass of the observable white dwarf in WD 0957- 



666 sy stem al lows it to be a hybrid white dwarf ([ben & 



Tutukc v 1985 ). Such white dwarfs descend from stars with 
initial mass 2.3 Mq - 5 Mq, which fill their Roche lobes 
in the stage of hydrogen burning in a shell around a non- 
degenerate helium core, become hot subdwarfs in the core 
helium-burning stage, but don't experience envelope ex- 
pansion after the formation of a degenerate carbon-oxygen 
core. Their masses are between 0.33 and 0.8 Mq. The for- 
mation of hybrid white dwarfs was considered in the study 
of the population of white dwarfs by |Tutukov fc Yungclsoii 



(1993Qiee their scenario 3) and all their following studies 
In a scenario shown in Fig. |^ (right) we start with 
a system of 2.4 Mq and 2.0 Mq in a relatively close or- 
bit (ao w 37 -Rq). At the instant of Roche lobe overflow 
the primary has a deep convective envelope and we ap- 
ply Eq. @ to compute the change in the orbital period. 
The primary becomes a compact helium star. A pecu- 
liarity of low-mass helium stars is their long life time, 



1.1 x 10 7 A/~ 3 - 75 yr for 0.33 < M/M < 0.7 ( |Pols ct al 
199l|)7|comparable to the lifetime of their main-sequence 
progenitors. As a consequence, the initially slightly less 
massive secondary fills its Roche lobe and becomes a he- 
lium white dwarf while the former primary still burns he- 
lium in its core. For some 250 Myr the system could be 
observed as a hot subdwarf with a companion unseen due 
to the difference in luminosities. After core helium exhaus- 
tion t he subdwarf cools a nd becomes a "hybrid" white 
dwarf. Moran et al. (1997 ) estimate from its T e g that the 



cooling age of this white dwarf is only ~ 10 7 yr and that 
the ratio of the luminosities of components is close to 5. 
This is compatible with the age of about 250 Myr expected 
in our scenario for the 0.31 Mq companion. 

7. Conclusion 

We followed the binary evolution for three double helium 
white dwarfs backwards and came to the following conclu- 



The last phase of mass transfer (the primary has al- 
ready become a white dwarf and the secondary fills its 
Roche lobe) was a spiral-in, for which we can constrain 
a, which describes the efficiency of orbital energy de- 
position into the common envelope, to lie between 
1 and 6, assuming a structure parameter A = 0.5. 
This efficiency value may be an overestimate since A 
may increase towards 1 at the end of the first red 
giant stage. Our result is in agreement with values 
of a ~ 4 found in population synthesis studies of 



low-mass X-ray binaries (Tauris 1996), double neutron 
stars (Portegies Zwart & Yungelson 1998) and double 
white dwarfs ( Nelemans et al. 2000] ). 
The parameters of all observed double helium white 
dwarfs may be reproduced with the same a A ~ 2. 
WD 0957-666 is the only system for which aA appears 
to be lower if both components are helium dwarfs. 
However, this system might have been formed with 
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aX w 2 if the immediate precursor of the currently ob- 
served white dwarf was a non-degenerate helium star 
and now it is a hybrid white dwarf. 

4. In order to explain the relatively high masses of the ob- 
served white dwarfs in close pairs, their direct progen- 
itors (i.e. white dwarf + giant binaries) must have had 
relatively wide orbits (between 25 and >500 days). The 
standard cases of mass transfer (Algol evolution and 
spiral-in) applied to the first phase of mass transfer, 
can not explain these intermediate wide orbits. Only 
if the masses of the observed white dwarfs are (much) 
lower than the current estimates (i.e. below 0.3 Mq) 
they could be formed through a phase of stable mass 
transfer. 

5. We suggest that in the first mass transfer phase for 
low-mass binaries with similar masses of the two stars, 
most of the mass of the envelope of the evolved star is 
lost without a significant spiral-in. This suggestion is 
supported by the fact that the original reasoning for 
spiral-in (drag forces in the envelope) is not applicable 
here, because the envelope can easily been spun up to 
corotate with the binary. 

In the absence of a physical picture for the removal 
of the envelope, we introduce a simple parameter 7 = 
(AJ/AMtot) (Mtot/J) to describe the loss of the an- 
gular momentum of the system as in the early work of 
jPaczyriski fc Ziolkowski (1967 ). Our analysis of the ob- 



served parameters for all observed double helium white 
dwarfs shows that the material of the envelope of the 
giant is expelled with 1.4 to 1.7 times the specific angu- 
lar momentum of the initial binary. The details of this 
kind of mass transfer should be investigated using 3D 
gas-dyna mical calculations, w hich are becoming avail- 
able (e.g. Bisikalo et al. 1998 ), but are not yet accurate 
enough to make predictions. 
6. With some well-constrained assumptions for the 
masses of the unseen companions in the other 8 double 
white dwarfs, we find similar results as for the double 
white dwarfs with two helium components. 
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